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High redshift Fermi blazars observed by GROND and Swift 
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ABSTRACT 

We observed 5 7-ray loud blazars at redshift greater than 2 with the X-Ray Telescope 
(XRT) and the Ultra Violet and Optical Telescope (UVOT) onboard the Swift satellite, and 
the Gamma-Ray burst Optical Near-Infrared Detector (GROND) instrument. These observa- 
tions were quasi simultaneous, usually within a few hours. For 4 of these blazars the near-IR 
to UV data show the presence of an accretion disc, and we could reliably estimate its accretion 
rate and black hole mass. One of them, PKS 1348+007, was found in an extraordinarily high 
IR-optical state, almost two orders of magnitude brighter than at the epoch of the Sloan Dig- 
ital Sky Survey observations. For all the 5 quasars the physical parameters of the jet emitting 
zone, derived by applying a one-zone emission model, are similar to that found for the bulk 
of other 7-ray loud quasars. With our observations we have X-ray data for the full sample 
of blazars at z > 2 present in the Fermi 2-yrs (2LAC) catalog. This allows to have a rather 
complete view of the spectral energy distribution of all high-redshift Fermi blazars, and to 
draw some conclusions about their properties, and especially about the relation between the 
accretion rate and the jet power 

Key words: galaxies: active-galaxies: jets-galaxies — radiation mechanisms: non-thermal 



^ ; 1 INTRODUCTION 



High redshift blazars are the most powerful persistent sources, and 
are usually connected with the most massive black holes (Ghis- 
ellini et al. 2009; 2010a; Volonteri et al. 2011). The Large Area 
Telescope (LAT) onboard the Fermi satellite (Atwood et al. 2009) 
and the Burst Alert Telescope (BAT) onboard Swift (Gehrels et al. 
2004) have provided tens of detections of blazars at 2: > 2 (Abdo 
et al. 2010, Ackermann et al. 2011; Aiello et al. 2009, Cusumano 
et al. 2010; Baumgartner et al. 201(o. The combinations of the 
two sets of data (Swift+Fermi) have allowed to measure the proper- 
ties and the bolometric luminosity of the jet non-thermal emission 
and to characterize the thermal component emitted by the accretion 
disc, namely the black hole mass Al and the accretion rate M. It is 
found (Ghisellini et al. 2010b; 201 1) that in powerful blazars black 
hole masses are usually greater than 10® M©, with disc luminosi- 
ties Ld '--^O.ll/Edd, jet kinetic powers Pj ~ Mc^, with Pj tightly 
related to Ld - This can test jet production models, since Pj ^ Mc^ 
requires that we are using another source of energy besides accre- 
tion, namely we must extract the black hole spin energy (see e.g. 
Tchekhovskoy, Narayan & McKinney 201 1). 
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In order to find the physical parameters of these sources it is 
very important to have a good coverage of their spectral energy dis- 
tribution (SED), that allows us to constrain the model parameters. 
Moreover, blazars are varying fast (hours to days, especially at high 
frequencies), and with large amplitudes (even by factor 10 or more 
in the 7-ray range, but sometimes even in the optical) implying that 
simultaneous observations are needed to well constrain their SED. 

High redshift blazars are interesting "per se", since we would 
like to know if and how the jet properties change with cosmic time. 
High redshift also usually means larger luminosities and powers, 
and this allows to study the more powerful jets. If we define a blazar 
as a source whose jet is observed with a viewing angle 9^ < l/F 
(F is the bulk Lorentz factor), then there must be other 2F^ similar 
objects whose jet is pointing elsewhere and whose flux is dramati- 
cally fainter (because of beaming). These misaligned sources share 
all the intrinsic properties of the blazar that is pointing at us, includ- 
ing the black hole mass. If we are able to estimate it for a blazar, we 
can put very interesting constraints on the density of heavy black 
holes of all radio sources in the young Universe (see Volonteri et 
al. 2011). 

From the z > 2 blazars detected by FennifLAT, we have cho- 
sen those with no data in the X-ray range (or with only an upper 
limit), and have organized a simultaneous observational campaign 
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Name 



RA 



Dec 



^7,1 



^7,2 



PKS 0519+01 
TXS 1149-084 
MG2 J133305+2725 
PMNJ1344-1723 
PKS 1348+007 



05 22 17.5 
11 52 17.2 
13 33 07.5 
13 44 14.4 
13 51 04.4 



+01 13 31 
-08 41 03 
+27 25 18 
-17 23 40 
+00 31 19 



2.941 

2.367 

2.126 

2.49 

2.084 



47.7 



48.5 



47.0 
47.7 
47.3 
48.2 
47.6 



Table 1. List of our sources. 1/7,1 refers to the [0.1-10 GeV] 7-ray lumi- 
nosity in the ILAC catalog, while 1/7,2 is the one in the 2LAC catalog, in 
units of erg s~^. 



involving the X-Ray Telescope (XRT) and the Ultra Violet and Op- 
tical Telescope (UVOT) onboard Swift (Gehrels et al. 2004) and 
the Gamma-Ray burst Optical Near-Infrared Detector (GROND) 
instrument (Greiner et al. 2008). The scientific rationale for observ- 
ing these blazars at X-ray and IR-optical frequencies is the follow- 
ing. In powerful blazars, the [0.3-10 keV] Swift band is where we 
expect the contribution of the inverse Compton emission of the jet. 
This is usually made by two components, according if the scattering 
process makes use of internally produced synchrotron seed photons 
(Synchrotron Self Compton, or SSC for short) or if the seeds are 
produced externally to the jet (External Compton, EC). They are 
usually characterized by a different spectrum and variability behav- 
ior, and the relative importance of the two gives information on the 
magnetic field and the bulk Lorentz factor. In the near IR, optical 
and UV bands, instead, we have the contributions of the accretion 
disc and the beamed synchrotron component. If we can distinguish 
the two contributions then we can estimate the black hole mass and 
the accretion rate (from the disc emission), and have information of 
the value of the magnetic field (from the synchrotron flux). 

The starting samples were the FennifLAT detected blazars at 
z > 2 present in the Fermi blazar catalogs after 1 1 months of op- 
erations (ILAC, Abdo et al. 2010) and after 2 years (2LAC, Acker- 
maim et al. 201 1). In the ILAC catalog there are 2 sources observed 
by Swift in 2009 for just a ks, for which we could only derive an 
upper limit to the X-ray flux. In the 2LAC catalog there are 3 new 
blazars at z > 2 (out of 31) without a proper characterization of 
their SED, because of no information on their X-ray flux, useful 
to characterize the jet beamed emission, nor good data coverage in 
the IR-optical-UV band, useful to derive the thermal (i.e. accretion 
disc) contribution. These 5 blazars are all at declination <+30°, and 
are visible from La Silla, where the GROND instrument operates. 
Therefore we organized quasi-simultaneous Swift and GROND ob- 
serving campaigns for these blazars. 

These observations provided us with optical-UV-X-ray in- 
formation on the entire high redshift Fermi sample ("clean" 2LAC 
with z > 2). Tab.[T]lists the 5 selected blazars, together with their 
7-ray k-corrected luminosities. These are the averaged luminosity 
over 11 months (ILAC) and 2 yrs (2LAC). For the blazars 1149- 
084 and 1344-1723, present in both catalogs, we give the corre- 
sponding two luminosities. 

With a good simultaneous coverage from the near IR to X-ray 
range, in addition to the Fermi data, we can properly characterize 
the SED in order to disentangle the non-thermal jet contribution 
and the thermal component, to find the black hole mass, the ac- 
cretion rate, and the physical jet quantities (magnetic field, bulk 
Lorentz factor, particle densities and jet power). 

We use a flat cosmology with Ho — 70 km s^^ Mpc^^, 
r2M=0.3 and the notation Q — 10^ Qx in cgs units. 



2 GROND OBSERVATIONS AND DATA ANALYSIS 

The 7-band GROND imager is mounted at the 2.2 m MPG/ESO 
telescope at La Silla Observatory (Chile). GROND is able to ob- 
serve simultaneously in 7 filters, from the NIR Ks (2300 nm) to 
the g' band (360 nm). Therefore it nicely complements the UVOT 
filters, with the bluest (g' and r') filter overlapping in part with the 
reddest v and b UVOT filter (and this is useful for cross calibration). 

We carried out observations for all sources simultaneously in 
all 1 g',r',i',z',J, H, Ks bands. The log of the GROND obser- 
vations and the related observing conditions are reported in Tab. 

m 

The GROND optical and NIR image reduction and photome- 
try were performed using standard IRAF tasks (Tody 1993), similar 
to the procedure described in Kruhler et al. (2008). A general model 
for the point-spread function (PSF) of each image was constructed 
using bright field stars, and it was then fitted to the point source. 
When the source field was covered by the SDSS (Smith et al. 2002) 
survey (i.e. PKS 0519+01, PKS 1348+007, MG2 J133305+2725), 
the absolute calibration of the g' ,r' ,i' , z bands was obtained with 
respect to the magnitudes of SDSS stars within the blazar field. In 
the other cases (i.e. PMN J1344-1723 and TXS 1149-084), opti- 
cal photometric calibration was performed relative to the magni- 
tudes of six secondary standards in the blazar field. During pho- 
tometric conditions, a primary SDSS standard field was observed 
within minutes of an observation of the source field. The obtained 
zero-points were corrected for atmospheric extinction and used to 
calibrate stars in the blazar field. The apparent magnitudes of the 
sources were measured with respect to these secondary standards. 
For all sources the J, A's bands calibrations were obtained with 
respect to magnitudes of the Two Micron All Sky Survey (2MASS) 
stars (Skrutskie et al. 2006). 

Tab.[3]reports the observed AB magnitudes, not corrected for 
the Galactic extinction listed in the last column and taken from 
Schlafly & Finkbeiner (201 1). 



3 SWIFT OBSERVATIONS AND DATA ANALYSIS 

We have analysed the Swift X-Ray Telescope (XRT; Burrows et 
al. 2005) and Optical-Ultraviolet Telescope (UVOT; Roming et al. 
2005) data. The data were screened, cleaned and analysed with the 
software package HEASOFT v. 6.12, with the calibration database 
updated to 22 March 2012. The XRT data were processed with the 
standard procedures (XRTPIPELINE v. . 12 . 6). Afl sources 
were observed in photon counting (PC) mode and grade 0-12 (sin- 
gle to quadruple pixel) were selected. The channels with energies 
below 0.3 keV and above 10 keV were excluded from the fit and the 
spectra were rebinned in energy so to have at least 20-30 counts per 
bin in order to apply the test. When there are no sufficient counts 
we applied the likelihood statistic in the form reported by Cash 
(1979). Each spectrum was analysed in XSPEC v. 12.7.1 with an 
absorbed power law model with a fixed Galactic column density as 
measured by Kalberla et al. (2005). The computed errors represent 
the 90% confidence interval on the spectral parameters. Tab. |4]re- 
ports the log of the observations and the best fit results of the X-ray 
data with a simple power law model. The X-ray spectra displayed 
in the SED have been rebinned to ensure the best visualization. 

UVOT source counts were extracted from a circular region 
5"-sized centred on the source position, while the background was 
extracted from an annulus with internal radius of 7" and variable 
outer radius depending on the nearest contaminating source. Data 
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were integrated with the uvotimsum task and then analysed by 
using the uvotsource task. The observed magnitudes have been 
dereddened according to the formulae by Cardelli et al. (1989) and 
converted into fluxes by using standard formulae and zero points. 

As can be seen, the UVOT observations yielded mostly up- 
per limits, and very few detections. The listed UVOT and GROND 
magnitudes do not take into account any difference between the 
two instruments, that is however likely, and of the order of 0.1- 
0.3 magnitudes (see Rau et al. 2012), especially if the observations 
were not exactly simultaneous, but separated by a few hours (the 
maximum separation - two days - occurs for PKS 0519+01). 



4 SPECTRAL ENERGY DISTRIBUTIONS AND 
MODELLING 

In Fig.[T}{5]we show the spectral energy distribution (SED) of our 
sources. We complement our near-simultaneous data with archival 
data taken from NED and ASDCB We show the Fermi/LAT data of 
the ILAC or 2LAC catalogs, but also the average flux correspond- 
ing to one month of Fermi observation, starting 2 weeks before and 
ending two weeks after the GR0ND-F5wi/f observing time. 

The adopted model is a one-zone and leptonic model, fully 
described in Ghisellini & Tavecchio (2009). The main properties 
are summarized in the Appendix, mainly to explain the meaning of 
the parameters listed in Tab.[6]and Tab. [7] 

In brief, the model assumes that the bulk of the jet dissipation 
takes place in one zone located at some distance i?diss from the 
black hole. For simplicity, the emitting region is assumed spherical 
with a radius R = 7/;_Rdiss, with = 0.1. The region is moving 
with a bulk Lorentz factor F, and is observed under a viewing angle 
0v Energetic electrons are injected throughout the source for a time 
equal to the light crossing time R/c, and the particle distribution is 
calculated (through the continuity equation) at this time, consid- 
ering radiative losses and possible electron-positron pair produc- 
tion and their reprocessing. In the following we briefly discuss the 
guidelines for the choice of the main parameters needed for the 
model. 



4.1 Guidelines for the choice of the parameters 

The luminosity of the accretion disc — It can be estimated di- 
rectly if the disc is visible and its spectrum peaks in the observed 
frequency range. This occurs for PKS 0519+01, for TXS 1149- 
084 and most likely for PMN 11344-1723, even if in this latter 
blazar there is a strong "contaminating" synchrotron component. 
The overall luminosity Ld of a standard accretion disc is roughly 
twice its uLiy) peak, therefore we can directly estimate Ld if we 
see a thermal peak in the SED. For MG2 J1333+2725 the IR to UV 
continuum is dominated by the steep tail of the synchrotron flux, 
while in PKS 1348+007 we have a hint of the contribution from 
the accretion disc from the photometric data of the SDSS. If an op- 
tical spectrum is available (as in the case of TXS 1 149-084 and 
PMN J1344-1723, Shaw et al. 2012), we have additional infor- 
mation from the luminosity of the broad emission lines. Through 
the templates of Francis et al. (1991) and/or of Vanden Berk et 
al (2001) we can reconstruct, from the luminosity of one or more 
lines, the entire luminosity of the broad line region, Lblr . Then, 



applying a typical covering factor C (namely, C ~ 0.1), we can 
estimate Ld- Finally, if no spectrum is available and the thermal 
component is completely swamped by the synchrotron flux, we can 
have a (rough) indication of Lblr through the correlation between 
Lblr and L-y as found by Sbarrato et al. (2012a). This has the form 

0.93 



L 



4L" 



(1) 



Note, however, that since can vary even by two orders of mag- 
nitude in the same source, the correlation necessarily has a large 
scatter, making the estimate of Lblr (and thus of Ld) uncertain. 

Black hole mass — If Ld is determined reliably, there is only one 
black hole mass value that can fit the flux produced by (a standard) 
accretion disc. In this case the derived mass is robust, with an as- 
sociated uncertainty of less than a factor 2 (see Fig. 1 in Sbarrato 
et al. 2012b, §4.4 and Fig.[T]below). If Ld is uncertain, this reflects 
also on the uncertainty on the derived black hole mass. 

Location of the emitting region — One of the specific features of 
the model is that it calculates the energy densities (magnetic and 
radiative) as a function of the distance i?diss from the black hole. In 
particular, if Rdiss < i?BLR, the energy density of the line photons 
as seen in the comoving frame becomes (up to a factor of order 
unity): 

l2lT 



Le 



AuRl 



(2) 
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where we have used -Rblr ~ '^-^^45 '-n*' ^ ~ 0.1. A 
similar relation holds for the IR radiation reprocessed by the torus 
located at a distance -Rir, namely when i?BLR < i?diss < -Rir- 
This limits 7?diss- 

Magnetic field — The magnetic energy density, although is for- 
mally a free parameter, must satisfy the Compton to synchrotron 
luminosity ratio, i.e. Lc/Lsyn = U^^d/U'B- t^rad includes both 
internally produced radiation (i.e. by synchrotron) and radiation 
produced externally (directly by the disc or reprocessed and re- 
isotropized by the BLR and the torus). 

Bulk Lorentz factor — The value of F determines the value of the 
radiation energy density of the external seed photons (oc F^) and 
hence the value of the magnetic field required to have the observed 
synchrotron to inverse Compton luminosity ratio. Further informa- 
tion come from the peak frequencies of the synchrotron and inverse 
Compton components, that depend also on the break energy of the 
electron distribution. 

Injected power — The power is injected throughout the source in 
the form of relativistic electrons. Through the continuity equation 
we calculate the particle distribution as a result of injection, cooling 
and possible pair production. The total injected power is such that 
the radiation produced by these particles agrees with the observed 
data. The injected distribution is assumed to be a smoothly broken 
power law (see Eq. |3]in the Appendix). The resulting distribution, 
modified by cooling, must agree with the observed slopes. 

4.2 Caveats 

Within the framework of the adopted model, there is some degen- 
eracy between a few set of parameters, that can be broken if some 
additional information, besides the SED, is provided. For instance, 
the bulk Lorentz factor and the viewing angle, together with the in- 
jected power in relativistic particles, can have a range of possible 



4 Ghisellini et al. 



Name Date Start time Exp: opt/IR Average seeing Average airmass 





yyyy-mm-dd 


[UTC] 


[s] 


[arcsec] 




0519+01 


2012-02-26 


00:14:00 


919/960 


0.70 


1.18 


1 149-084 


2012-03-20 


02:09:51 


1501/1200 


0.72 


1.27 


1333+2725 


2012-05-09 


02:40:11 


426/720 


1.31 


1.83 


1344-1723 


2012-04-04 


08:36:13 


919/960 


1.32 


1.44 


1348+007 


2012-03-16 


05:03:37 


3002/2400 


1.15 


1.26 


1348+007 


2012-03-18 


09:38:11 


460/480 


1.43 


1.61 



Table 2. Log of the GROND observations. Exposures refer to optical/NIR filters while the average seeing is calculated in the r' band. 



Name 


9' 


r' 


i' 


z' 


J 




H 




Av 


0519+01 


20.37±0.05 


19.83±0.05 


19.50±0.05 


19.28±0.05 


18 


..74±0.11 


18.54±0.11 


18.6I±0.15 


0.38 


1149-084 


19.65±0.05 


19.42±0.06 


19.40±0.08 


19.04±0.08 


18 


..97±0.11 


18.67±0.12 


18.24±0.13 


0.23 


1333+2725 


20.38±0.06 


19.89±0.06 


19.39±0.06 


19.07±0.05 


18 


..33±0.10 


17.73±0.10 


17.44±0.13 


0.03 


1344-1723 


20.69±0.05 


20.26±0.05 


19.99±0.05 


19.56±0.06 


19.01±0.12 


18.42±0.12 


18.11±0.13 


0.37 


1348+007 (1) 


19.65±0.05 


19.25±0.05 


18.90±0.05 


18.62±0.06 


18 


..25±0.10 


17.79±0.10 


17.40±0.12 


0.11 


1348+007 (2) 


20.24±0.05 


19.88±0.05 


19.49±0.05 


19.19±0.06 


18 


.66±0.11 


18.20±0.11 


17.76±0.12 


0.11 



Table 3. Observed magnitudes, not corrected for Galactic foreground reddening, in the AB system. Errors include systematics. The last column reports the 
value of the Galactic Ay from Schlafly & Finkbeiner (2011). 



values. If, in addition, we have a limit on the variability timescale 
and/or the superluminal speed, then we can chose a unique set of 
parameters. 

The black hole mass found by fitting a Shakura-Sunjaev disc 
to the near IR-optical-UV data gives excellent results if the flux 
at these frequencies is not contaminated by the synchrotron flux. 
Otherwise, the disc luminosity can be estimated rather accurately 
from the observed broad line luminosities, but the black hole mass 
can have a large uncertainty, partly mitigated by assuming that the 
disc cannot be super-Eddington (thus yielding a lower limit to the 
black hole mass). 

Also the derived jet powers bear some uncertainties due to sev- 
eral unknowns: i) we do not know if charge neutrality is provided 
by protons, or by positrons. Recent studies (Ghisellini & Tavecchio 
2010) have shown that a pure pair plasma would suffer a severe 
Compton drag while crossing the broad line region, limiting the 
positron/proton ratio to nearly 20 (in agreement with independent 
estimates put forward by Sikora & Madejski 2000); ii) we can es- 
timate the amount of emitting particles, but there can be additional 
particles that are not accelerated, but nevertheless participate to the 
bulk motion of the jet, hence to its kinetic power; iii) the estimated 
magnetic field B is the one in the emitting region. If the dissipation 
mechanism is magnetic reconnection, it is likely that in the emitting 
region B is smaller than in the surroundings. 

Bearing in mind these limitations, we now discuss the physical 
parameters found by adopting our model. 

4.3 Physical parameters 

Tab.[6]lists the parameters of the applied model. It also reports the 
parameters adopted in Ghisellini et al. (2011) for the two sources 
in common with that paper. Consider that the SED available at the 
time of Ghisellini et al. (2011) was largely incomplete, lacking the 
WISE, GROND and Swift data. Moreover, also the high frequency 
radio coverage has been improved with data provided by the Planck 



and WMAP satellites (see also Giommi et al. 2012 for a collection 
of blazar's SED including Planck and WMAP data). The improved 
characterization of the SED of these blazars allowed a better esti- 
mate of the physical parameters: we find that although the values 
found now are not very different from what we have guessed be- 
fore, the uncertainty is much less. 

In general, all the derived parameters are well within the distri- 
butions derived for a large sample of 7-ray loud blazars studied in 
Ghisellini et al. (2010a). For all 5 blazars in our sample, the region 
dissipating most of the flux we see is located at several hundreds 
of Schwarzschild radii from the black hole, with bulk Lorentz fac- 
tors in the range 10-15, and small viewing angle {6^ ~ 3°). The 
magnetic field is in the range 1-8 Gauss, and the intrinsic power 
injected in the form of relativistic electrons is of the order of 10^'* 
erg as measured in the comoving frame. The black hole mass 
is around M ~ 10® M©, with a luminosity of the accretion disc 
ranging from lO^'^ erg s"^ (for PKS 1348+007) to ~ (2-3) x 
erg s" ^ (for PKS 05 19+01 and TXS 1 149-084). These are all very 
typical values for FSRQs. In the following we discuss in more de- 
tail individual sources. 

4.4 PKS 0519+01 

Present in the 2LAC catalog, this 7-ray blazars had no previous X- 
ray information. There is no redshift in NED, the value for z comes 
from the value listed in the 2LAC catalog, but the spectrum is still 
unpublished. Our observations (and the WISE data) greatly improve 
our knowledge of the SED, as shown in Fig.[T] despite the fact that 
for this source, the UVOT data were only upper limits. Note that for 
each UVOT upper limit we plot two arrows, corresponding to the 
observed datum and the de-absorbed one. The latter is derived as- 
suming, along the line of sight, an average distribution of absorbing 
Lya clouds (see Fig. 3 in Ghisellini et al. 2010a). 

The GROND photometric points allow to determine the peak 
of the accretion disc spectrum, hence its luminosity Ld ^ 2x Iff"^ 
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Name 


Date 




Exp. 


luGal 


r 


Fx 


(d.o.f.) 


Note 




yyyy-mm-ddAJTC) 




ks 






10-1^ cgs 






0519+01 


2012-02-24/01:26 




21.9 


1.07e21 


1.7640.24 


3.1 


76.28 (95) 


Cash Stat 


1 149-084 


2012-03-19/00:10 + 


2012-03-20/08:09 


17.1 


4.75e20 


1.4940.28 


2.5 


70.02 (67) 


Cash Stat 


1333+2725 


2012-05-09/00:09 4 


-2012-05-10/21:16 


19.8 


1.02e20 


1.8440.17 


6.65 


10.87(12) 




1344-1723 


2012-04-03/10:55 4 


-2012-04-04/01:41 


14.9 


8.70e20 


1.8740.40 


1.5 


36.68 (39) 


Cash Stat 


1348+007 


2012-03-16/04:51 4 


-2012-03-18/13:17 


17.5 


2.22e20 


1.5040.26 


2.9 


82.27 (78) 


Cash Stat 



Table 4. Observation log and results of the X-ray analysis of XRT data. The flux Fx is in the [0.3-10 keV] band, it is de-absorbed with the indicated N^^^ 
and it is in units of IQ-^"' erg cm~^ s~^. Apart from 0519+01, all other sources were observed in two occasions: the listed exposure time is the sum. The 
analysis has been performed on the total. All sources but 1333+2725 have been analyzed with the Cash statistics (Cash 1979; Gehrels 1986). 



Name 


V 


b 


u 


uvwl 


uvm2 


uvw2 


0519+01 


>19.1 


>20.0 


>19.5 


>19.6 


>21.7 


>21.9 


1 149-084 


>19.1 


>20.0 


>19.7 


>20.0 


>20.6 


>21.8 


1333+2725 


>I9.2 


>20.2 


20.3 40.3 


20.440.4 


19.740.2 


19.940.1 


1344-1723 


>19.0 


>20.0 


21.040.2 


>19.9 


>20.1 


>20.7 


1348+007 


>19.1 


>20.2 


>20.8 


>21.0 


20.3 40.1 


>20.8 



Table 5. UVOT magnitudes. The magnitude lower limits are at the 3cr level. 



erg and black hole mass, that turns out to be M — 4.5 x 
10® M0. This is the largest value we find for the 5 blazars here 
considered. To estimate the uncertainty for this value, we show in 
Fig. [lithe fit with a black hole mass of 2.25, 4.5 and 9 billion solar 
masses, as labelled, keeping the same id. The fit with the largest 
mass underestimates the high frequency GROND point, while the 
fit with the lowest mass underestimates all but the high frequency 
GROND fluxes. We can conclude that the mass determination has 
an uncertainty, in this case, of less than a factor 2. 

The bottom panel of Fig. [T] shows the complete SED of the 
source together with the model. We show, separately, the syn- 
chrotron component (solid light green line) the torus+disc+X- 
ray corona contribution (black short-dashed line) and the inverse 
Compton contribution, dominated by scattering with emission line 
seed photons (grey dot-dashed line). The thick (blue) solid line is 
the sum. We show also the Fermi upper limit on the 7-ray flux re- 
sulting from a month of data centered on the time of the Swift and 
GROND observations (thick red arrow). 



4.5 TXS 1149-084 

This source has been observed spectroscopically in the optical by 
Shaw et al. (2012), that reported a luminosity (Lciv = 2 x 10'*'' 
erg s"^) and a FWHM (7200 km s"^) of the CIV broad emis- 
sion line, together with the luminosity of the continuum at 1350 
A (ii35o = 1.3 X lO^'^erg s"^). These data allowed Shaw et al. 
(2012) to estimate a black hole mass of M = 2.4 x 10^ Mq, ap- 
plying the virial method. Using the template of Francis et al. (1991) 
and Vanden Berk (2001), one can derive the overall luminosity of 
the BLR, Lblr ~ 2 x 10*^ erg s~^ and then a disc luminosity 
ten times greater (assuming a covering factor equal to 0.1). This 
value agrees well with the GROND+Swi/? data, from which we de- 
termined the peak of the disc component, with Ld = 3.2 x 10''® 
erg s~'. Therefore also in this case the black hole mass is well de- 
termined, A/ = 1.5 X 10® M0. 



In this source the WISE data, together with the high frequency 
radio data (from Planck) show a strong synchrotron component, 
peaking in the submm range. This is also indicated by the GROND 
data, showing un upturn towards the low frequencies. This upturn 
constrains the possible models capable of reproducing the syn- 
chrotron peak. The self-absorption frequency of our compact emit- 
ting zone occurs at ~ 760 GHz (observed frame), making the syn- 
chrotron component very narrow. We derive a rather large magnetic 
field (~8 G), to account for the strength of the synchrotron flux. Ac- 
cordingly, also the synchrotron Self Compton flux is not negligible, 
and contributes to the soft X-rays (long dashed grey curve in Fig. 

13. 

Also for this source we show the FermifLAT upper limit mea- 
sured from 1-month of data around the time of the Swift and 
GROND observations. 

Tab.[6]lists also the parameters used in Ghisellini et al. (201 1), 
for which no GROND data were available and there was only an 
upper limit to the X-ray flux (shown as a blue arrow in Fig. |2]l, 
implying that the source has brightened in the X-ray band. The 
main differences with those results concern the black hole mass 
(it was ^3 times greater), the value of Rdiss (4 times larger) and 
the magnetic field (5.5 times smaller). The previous UVOT fluxes 
were slightly larger, and in the absence of additional optical IR data 
these resulted in an overestimation of the black hole mass and and 
disc luminosity (see Tab. |6]l, instead of a larger synchrotron flux. 
This well illustrates the importance of having a good coverage in 
the IR-optical, and also some information on the emission lines 
(from Shaw et al. 2012). The better coverage in the near and far 
IR and in the submm range allows to characterize better the syn- 
chrotron component, while the detection in the X-rays (at a flux 
larger than the previous upper limit) allows to determine the impor- 
tance of the inverse Compton process (sum of SSC and EC). We 
find that the overall SED can be explained assuming a relatively 
strong synchrotron (and SSC) components, consequence of a mag- 
netic field larger than the one assumed in Ghisellini et al. (201 1). 
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Name 


z 


-^diss 


M 








B 


r 


7b 






S2 


[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


[7] 


[8] 


[9] 


[10] 


[11] 


[12] 


[13] 


0519+01 


2.941 


540 (400) 


4.5e9 


486 


0.011 


23.6 (0.035) 


1.9 


11 


50 


3e3 


1 


2.5 


1 149-084 


2.367 


180(400) 


1.5e9 


561 


0.01 


31.5 (0.14) 


7.9 


15 


300 


800 


-1 


3.3 






720 (600) 


4e9 


849 


0.015 


72(0.12) 


1.39 


14 


300 


3e3 


-1 


2 


1333+2725 


2.126 


45 (1500) 


le8 


300 


0.01 


9 (0.6) 


5.5 


13 


600 


1.8e3 





2.4 


1344-1723 


2.49 


225 (500) 


1.5e9 


351 


0.021 


12.3 (0.055) 


2.3 


13 


1.4e3 


2.3e3 





2.7 






330 (1100) 


le9 


274 


0.027 


7.5(0.05) 


0.89 


13 


1.4e3 


Se3 


-1 


2.5 


1348 high 


2.084 


60 (500) 


4e8 


115 


8.5e-3 


1.3(0.022) 


4.5 


11 


500 


4.7e3 


-1 


2.2 


1348 low 


2.084 


96 (800) 


4e8 


115 


2.2e-3 


1.3 (0.022) 


1.3 


15 


300 


3e3 


-1 


2.5 



Table 6. List of pai'ameters used to construct the theoretical SED. Not all of them are "input parameters" for the model: /?blr, is uniquely determined from 
Ld. and the cooling energy 7c and U' are derived parameters. Col. [1]: name; Col. [2]: redshift; Col. [3]: dissipation radius in units of 10^^ cm and (in 
parenthesis) in units of Schwarzschild radii; Col. [4]: black hole mass in solar masses; Col. [5]: size of the BLR in units of 10^''' cm; Col. [6]: power injected in 
the blob calculated in the comoving frame, in units of 10*^ erg s^^; Col. [7]: accretion disc luminosity in units of 10^^ erg s^^ and (in parenthesis) in units of 
^Edd; Col. [8]: magnetic field in Gauss; Col. [9]: bulk Lorentz factor at /?diss; Col. [10] and [1 1]: break and maximum random Lorentz factors of the injected 
electrons; Col. [12]: and [13]: slopes of the injected electron distribution [Q(7)] below and above 7]-,; For all sources we assumed a viewing angle 6v = 3°. 
The total X-ray corona luminosity is assumed to be in the range 10-30 per cent of L^j. Its spectral shape is assumed to be always oc exp(— /ii//150 keV). 
The parameters in italics refer to the physical quantities found in Ghisellini et al. (201 1). 



Name 


logPr 




logPc 


logPp 


logPj, 


0519+01 


45.00 


45.68 


44.53 


46.94 


45.30 


1 149-084 


45.08 


46.22 


43.93 


46.16 


45.38 




45.47 


45.87 


43.83 


46.24 


45.77 


1333+2725 


45.18 


44.58 


44.33 


45.96 


45.48 


1344-1723 


45.54 


45.23 


44.07 


46.12 


45.84 




45.66 


44.74 


44.43 


46.03 


45.96 


1348 high 


45.00 


44.52 


43.96 


45.38 


45.30 


1348 low 


44.68 


44.12 


43.59 


45.34 


44.98 



black hole mass we have no strong constraints but a slight flattening 
of the GROND data towards the blue, possibly indicating an upturn 
of the spectrum. If due to the presence of the accretion disc, this 
implies a relatively small mass (i.e. a high maximum temperature), 
so we have chosen an illustrative value of M — IO^A^q. For this 
mass the disc must emit at 60% the Eddington rate. 

A strong synchrotron component implies that the SSC flux can 
contribute to the soft X-ray spectrum, as shown in Fig. |3] while 
the EC component dominates the bolometric output as in the other 
sources. 



Table 7. Logarithm of the jet power in the form of radiation (Pi), Poynting 
flux (Pb), bulk motion of electrons (Pc) and protons (Pp, assuming one 
proton per emitting electron). The last column lists the minimum jet power, 
calculated assuming that the radiation drag of the jet halves its bulk Lorentz 
factor This limit corresponds to twice the radiated power Pr . Powers are in 
erg s~^. The parameters in italics refer to the physical quantities found in 
Ghisellini et al. (2011). 

4.6 MG2 J133305+2725 

This is a blazar that is present in the photometric optical SDSS 
survey (with a magnitude r — 20.18), but with no spectroscopic 
observations. UVOT detected the source in the bluest filters. At 
these frequencies the flux is partially absorbed by intervening Lya 
clouds, whose total optical depth can be roughly calculated by av- 
eraging over many line of sights (as done in Ghisellini et al. 2010). 
There is however a large dispersion around these mean values, and 
for this reason we have indicated, with bars, the possible range of 
the de-absorbed UVOT fluxes. 

The slope defined by the GROND data is steep, not consistent 
with a disc spectrum, that must therefore be hidden by the tail of 
the synchrotron flux. A strong synchrotron component is indeed 
needed by the WISE, WMAP and Planck data. However, for the 
fit, we have given less weight to the WISE data, since they are not 
simultaneous. 

The emission disc component shown in Fig.[3]is only illustra- 
tive, but its luminosity cannot much be less than shown, due to the 
presence of the broad emission lines that make this blazar a FSRQs 
(but there is no published information on the line strength). For the 



4.7 PIVIN J1344-1723 

This source has been observed spectroscopically in the optical by 
Shaw et al. (2012), that derived the luminosity of the CIV broad line 
(Lciv ~ lO*** erg s~^), corresponding (using Francis et al. 1991 
or Vander Berk et al. 2001), to I/blr ~ lO"*^ erg s"^ Adopting 
a covering factor C ~ 0.1 we then have Ld ~ 10*'' erg s^^. The 
GROND data show a flattening (in uFv) of the spectrum, that we 
interpret as the emergence of the accretion disc component. With 
UVOT we have upper limits in all filters except in the it band. 
Note that in u the source is already affected by possible absorp- 
tion by intervening Lya clouds. In Fig. |4] both the observed and 
the de-absorbed flux are plotted. However, we caution about the 
large uncertainty connected with the use of an average distribution 
of absorbers along the line of sight. By decomposing the optical- 
UV emission with a synchrotron+accretion disc component, we de- 
rived Ld ~ 1.2 X 10*^ erg s"'^, in agreement with what is in- 
dicated by the CIV broad line. We then derive a black hole mass 
Al — 1.5 X 10^ Mq. This value is about the same as that derived 
by Shaw et al. (2012) using the virial method, the FWHIVI of the 
CIV broad line (6000 km s"'^) and the continuum luminosity at 
1350 A (L1350 = 5.9 X lO^'^erg s'^), giving M = 1.3 x IO'^Mq. 

The WISE data indicate a strong synchrotron component. 
However, these data cannot connect smoothly with the GROND IR 
points, strongly suggesting that the synchrotron emission is vari- 
able with a large amplitude. This implies also that the blazar was in 
a somewhat low state at the epoch of our observations, and for this 
reason we did not attempt to accurately reproduce the 7-ray flux. 
For the model shown in Fig. |4]we have derived 2.3 Gauss for the 
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Figure 1. Top panel: a zoom of the SED of the blazar PKS 05 19+01 in the 
IR, optical and UV. Light (cyan) arrows are the upper limit as observed by 
UVOT, darker (red) arrows are the same data de-absorbed by the predicted 
amount of Lyo forest absorption along the line of sight (see Ghisellini et 
al. 2010a). Triangles (orange) are the GROND data, empty circles are the 
WISE data. We also show three accretion disc models (short dashed lines) 
with the same accretion luminosity and three different black hole mass, 
from 2.25 to 9 billion solar inasses (as labelled), together with the con- 
tribution of the torus emission, emitting in the IR. The thin solid (green) 
line is the synchrotron component, and the solid (cyan, black and red) lines 
are the sum of the accretion disc+torus+synchrotron flux. Note that WISE 
data are not simultaneous. Bottom: The entire SED of the PKS 0519+01, to- 
gether with the adopted model. Filled (red) circles are the XRT data, and the 
(heavy) arrow in the 7-ray band is the upper limit con'esponding to 1 month 
of Fermi data centered to the SwiTf+GROND observations, while the other 
arrows in the 7-ray band correspond to data from 2yr integration. Note also 
the detection at ~ 10^'' Hz. Archival data are form the online service of 
ASI Science Data Center (ASDC, green filled circles) and NED (empty cir- 
cles). Short dashed Une: contribution from the accretion disc (with a black 
hole mass M = 4.5 X 10^ Mq), IR torus and corona. Thin solid (green) 
line: synchrotron; dot-dashed line: EC component. Solid (blue) line: sum 
of all components. 



Figure 2. Lines and symbols as in Fig.[T] for TXS 1 149-084. Bottom panel: 
in this case the SSC component (long dashed grey line) contributes to the 
soft X-ray flux. The upper limit in X-rays (blue arrow) corresponds to an 
earlier and very short ( 1 ks) Swift observation, when the source was detected 
in two UVOT filters as labelled in the top panel ["UVOT (previous)"]. 

magnetic field: with this value the SSC barely contributes to the soft 
X-rays. At higher X-ray energies the flux is completely dominated 
by the external Compton process (with line photons as seeds). The 
present set of data can be compared with what was known previ- 
ously, and studied in Ghisellini et al. (2011, see their Fig. 5). One 
can appreciate the great improvement, and consequently the im- 
proved confidence on the derived physical quantities. 

4.8 PKS 1348+007 

This blazar showed an extraordinary optical flare, as derived by 
comparing the optical SDSS data (r = 22.4) with our GROND and 
UVOT data, together with the IR flux seen by WISE. Fast variability 
is present also in our GROND data, taken two nights apart (see Tab. 
|2l Tab.[3]and Fig.|5j. The source varied by about half a magnitude 
in all filters (except Ks), i.e. by ^^60% in flux. There is also a hint 
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Figure 3. Lines and symbols as in Fig.[T] for MG2 J133305+2725. Also 
in this case, as for TXS 1 149-084, the SSC component contributes to the 
soft X-ray flux. The bars in the bluest UVOT filters indicate the possible 
range of intrinsic flux levels: the lowest extremes correspond to the detected 
flux, and the highest extremes to the flux once de-absorbed by the average 
intervening Lyo absorption as calculated in Ghisellini et al. (2010). 



(although marginal) of a "harder when brighter" behaviour. Unfor- 
tunately, there are no SDSS spectra for this source. Comparing our 
data with the SDSS photometric data, obtained on May 20, 2009, 
the synchrotron emission had to vary by a factor ~30. As can be 
seen in Fig. [5] the synchrotron component completely outshines the 
disc emission at the time of our GR0ND+5wi// observations. For- 
tunately, the SDSS data hint to the presence of an accretion disc, 
through a flat (in vF^) slope suggested by the photometric data. 
We have assumed that the SDSS fluxes are completely produced 
by the accretion disc, and derived a luminosity Ld — 1.3 x 10''^ 
erg s^^ and a black hole mass M — A x 10* A/©. 

The extraordinary optical variability of PKS 1348+007 is not 
unprecedented, being similar to the optical flare shown by 3C 454.3 
in 2005 (Fuhrmann et al. 2006; Pian et al. 2006; Villata et al. 2006; 
Giommi et al., 2006). In principle, these large optical variations 



Figure 4. Lines and symbols as in Fig.[T] for PMN J1344-1723. 



could be explained in two different scenarios. In the first, one can 
assume that the jet power is unchanged, but the dissipation regions 
varies and is smaller for the high optical state. The magnetic field, 
which should decrease along the jet, is larger for smaller i?diss, 
implying a larger U'g . On the other hand, the radiation energy den- 
sity of the broad line photons is constant within i?BLR, and in the 
comoving frame is ?7blr ~ r^/(127r) (see Eq. 2). Therefore the 
synchrotron to inverse Compton luminosity ratio Lsy-n/Lc changes 
for varying i?diss even if the jet carries the same amount of power. 
In this case the 7-ray luminosity could remain constant or even de- 
crease during an optical flare. Alternatively, the jet power can vary, 
making the optical and the 7-ray fluxes vary together. In this case 
an high optical state should be accompanied by a larger 7-ray flux. 

Unfortunately, we do not have information about the high en- 
ergy emission for the low optical state of the source, so we can- 
not distinguish between these two hypotheses. We have simply 
explored the first option, looking for a solution for both states 
that maintains the total jet power roughly constant. The two mod- 
els shown in Fig. |5] correspond to the jet emission produced at 
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Figure 5. Lines and symbols as in Fig.[T] for PKS 1348+007. Note tlie large 
amplitude variability in the optical. The SDSS photometric fluxes, taken on 
May 20. 2009, are a factor ~30 below the GROND+Swi/i'+WISE data. Fast 
variability is present also in our GROND data, taken two nights apart (see 
Tab.[2] Tab.[3). 



^diss = 60 and 96 Schwarzschild radii, with F = 11 and 15, 
respectively, and with a equal power in bulk motion of the cold pro- 
tons. The magnetic field is 4.5 G for the high optical state (smaller 
i?diss and smaller F) and 1.3 G. for the low optical state (larger 
i?diss and larger F). The results demonstrate that this case is indeed 
possible. Simultaneous observations in the X-rays and in the 7-ray 
range when the source is in a low optical state can indeed decide if 
this is what really occurs. 



5 DISCUSSION AND CONCLUSIONS 

With our 5w;/f+GR0ND observational campaign, we have secured 
the X-ray coverage for all the blazars at z > 2 present on the 
"clean" 2LAC catalog. The use of simultaneous GROND and Swift 
observations were crucial to find the black hole mass and accretion 
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Figure 6. The minimum jet power Pj ^in (top panel) and the total jet power 
Pj (including one proton per emitting electron), as a function of the ac- 
cretion luminosity Lj. Our sources (black stars) are compared with other 
7-ray loud FSRQs studied previously (BAT z > 2: Ghisellini et al. 2010a; 
GGIO: Ghisellini et al. 2010b; Blue QSOs: Ghisellini et al. 2012; S 12: Sbar- 
rato et al. 2012a). The blazars studied in this paper lie on the bulk of the 
distribution. The minimum jet power is twice the power that the jet spends 
to produce the radiation we see (i.e. fj,min = 2Pr, see text and Ghisellini 
& Tavecchio 2010). 



rate for 3 out of 5 sources (for MG2 J 133305+2725 the synchrotron 
jet component was too strong to see the accretion disc emission, 
and for PKS 1348+007 our observations cought the source in a 
very high state, hiding the disc emission that was instead visible 
by previous SDSS photometric observations). We find that both M 
and Ld/^Edd are not extreme, but rather standard for the powerful 
FSRQs detected by FermifLAT. We have shown that the method 
of combining broad line luminosities, near IR/optical luminosi- 
ties (when the disc is visible), and a standard Shakura & Sunyaev 
(1973) disc emission model is very powerful to find M and M. 
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With good data, showing the disc emission at its peak, the un- 
certainty on the black hole mass is less than a factor 2. Since the 
blazars here studied are 7-ray emitters, we can also robustly con- 
strain the jet power, since the 7-ray luminosity, in these blazars, is 
almost equal to the bolometric one. We can then study in a robust 
way the link between the accretion and the jet powers. 

In §4.2 we discussed the uncertainties related to the jet power, 
associated to the unknown proton/lepton ratio. One robust lower 
limit is associated to Pr, the power that the jet spends to produce the 
radiation we see. It is simply Py ~ F'^Lboi (see Appendix), where 
Lboi is the total jet luminosity. Ghisellini & Tavecchio (2010) have 
discussed the importance of the Compton rocket effect on the jet 
when it is crossing the BLR. In the comoving frame of the jet, the 
seed photons are not isotropically distributed. This implies that also 
the inverse Compton scattered photons are not isotropic, and more 
power is emitted in the forward direction (i.e. along the jet veloc- 
ity direction). The jet then must recoil. If the jet is "heavy" (i.e. 
one proton per electron) the recoil is negligible, but if the jet is 
made by pairs the effect is very important. The jet halves its bulk 
Lorentz factor when there are ~20 pairs per proton. This corre- 
sponds to a power that is roughly equal to Pr. Requiring that the 
jet does not decelerate significantly, we end up with a minimum jet 
power (corresponding to a minimum amount of protons) that sim- 
ply is Pj,min = 2Pr. On the other hand, if we assume no pairs and 
therefore one proton per electron, we have a jet power Pj. 

These two quantities are listed in the two last columns of Tab. 
13 Fig.|6]shows Pj,miii (top panel) and Pj (bottom panel) as a func- 
tion of Ld for our blazars, where they are compared to other pow- 
erful FSRQs that we have analysed in the past. These are the FS- 
RQs at jz > 2 detected by the 3 years all sky survey of Swi/z/BAT 
(Ajello et al. 2009, these blazars are labelled BAT z > 2 in Fig.|6] 
and were studied by Ghisellini et al. 2010a); the FSRQs detected by 
FermifLAT in the first 3 months of operations (labelled GIO; Ghis- 
ellini et al. 2010b); the 4 "blue" quasars (FSRQs with a strong syn- 
chrotron component peaking in the optical, labelled "Blue QSOs" 
in Fig.m Ghisellini et al. 2012); and all the FSRQs in the ILAC 
sample present also in the SDSS spectroscopic survey (labelled 
S12, Sbarrato et al. 2012a). 

Note that: 

• Our blazars lie in the bulk of the distribution, with average 
values of the jet power and accretion luminosity. 

• The correlation between Pj,miii and Ld (top panel of Fig.|6l( is 
significantly less dispersed than the Pj-I/d relation. We re-iterate 
that Pj is found considering one proton per emitting electron, so 
that a non-constant number of pairs per proton could be responsible 
for the larger dispersion. However, we think it is premature to draw 
any strong conclusion, given the related uncertainties. 

• ^j.min is of the Same order as Ld- Given that this jet power is 
a lower limit, this suggest that the total jet power can be larger than 
Ld- In turn, this suggests that the origin of the jet power cannot 
be accretion only, and favours the extraction of the black hole spin 
energy as the prime movers of the jet power. 

• We can compare the z < 2 blazar detected by SwiftfBAT with 
the blazar in our sample. It is evident that the former have both 
more powerful jets and more luminous accretion discs, lying at the 
higher end of the distribution of powers. 

With our observational campaign all the Fermi z > 2 blazars 
in the 2LAC catalog have been observed and detected in the 
Swift/XKF energy range. This allows us to have a conclusive view 
of the SED of high redshift 7-ray blazars. 

From our results we conclude that Fermi blazars at high red- 



shifts are indeed powerful, but not extreme. Similarly, also their 
black hole masses are large, but not extreme. This can be con- 
trasted with high-2: blazars detected at hard X-ray energies, that 
all have extreme values of the jet power, of the disc luminosity and 
of the black hole mass. Therefore the hard X-ray band (>30 keV) 
is more efficient than the 7-ray band (given the current sensitivi- 
ties) in finding the most powerful blazars. This is expected if the so 
called "blazar sequence" (Fossati et al. 1998; Ghisellini et al. 1998) 
holds even at the highest power, since it predicts that the peak fre- 
quencies of both the synchrotron and the high energy humps shift 
to lower values when increasing the bolometric observed luminos- 
ity. At the highest end of the power distribution, the high energy 
peak can shift to sub-MeV energies, implying a large hard X-ray 
flux and a smaller 7-ray flux (with the K-correction working in the 
same direction). 
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APPENDIX 

At a distance Pdiss from the black hole of mass M the jet dissipates 
part of its power and injects relativistic electrons throughout the 
emitting region, assumed to be spherical, with radius R — i/; Pdiss, 
with tp — 0.1. In the region there is a tangled magnetic field B. The 
relativistic electrons are injected with a smoothly joining broken 
power law in energy: 



0(7) = Qo 



1 + (7/7b)-^i+^2 



(3) 



The energy particle distribution A'^(7) [cm^^] is calculated solving 
the continuity equation where particle injection, radiative cooling 
and pair production (via the 7-7 — )■ process), are taken into 
account. The created pairs contribute to the emission. 

The injection process lasts for a light crossing time P/c, and 
we calculate A''(7) at this time. This assumption comes from the 
fact that even if injection lasted longer, adiabatic losses caused by 
the expansion of the source (which is traveling while emitting) and 
the corresponding decrease of the magnetic field would make the 
observed flux to decrease. Therefore the calculated spectra corre- 
spond to the maximum of a flaring episode. 

The total power injected into the source in the form of rel- 
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ativistic electrons is P/ — nicC^V J (3(7)7^7, whiere V — 
(4tv/3)R^ is the volume of the emitting region. 

The bolometric luminosity of the accretion disc is Ld ■ Above 
and below the accretion disc, in its inner parts, there is an X-ray 
emitting corona of luminosity Lx (it is fixed at a level of 30% 
of Ld)- Its spectrum is a power law of energy index ax ~ 1 
ending with a exponential cut at Ec =150 keV. The specific en- 
ergy density (i.e. as a function of frequency) of the disc and the 
corona are calculated in the comoving frame of the emitting blob, 
and used to properly calculate the resulting External inverse Comp- 
ton spectrum. The BLR is assumed to be a thin spherical shell, of 
radius -Rblr ~ lO^^^^fs cm. We consider also the presence of a 
IR torus, at larger distances. The internally produced synchrotron 
emission is used to calculate the synchrotron self Compton (SSC) 
flux. Table[6]lists the adopted parameters. 

The power carried by the jet can be in the form of radiation 
(Pr), magnetic field (Pb), emitting electrons (Pc, no cold electron 
component is assumed) and cold protons (Pp, assuming one proton 
per emitting electron). All the powers are calculated as 



P = nR^r^lScUi 



(4) 



where f/j' is the energy density of the i component, as measured in 
the comoving frame. 

The power carried in the form of the produced radiation, Pr = 
nR^r^fic C/rad> can be re-written as [using U^g_d = L' / (47rP'^c)]: 

Pr = L'— = L— ~ L— (5) 
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where L is the total observed non-thermal luminosity (L' is in the 
comoving frame) and U^^^ the radiation energy density produced 
by the jet (i.e. excluding the external components). The last equality 
assumes 6^ ~ l/F. 

When calculating Pc (the jet power in bulk motion of emitting 
electrons) we include their average energy, i.e. = nc{'))maC? ■ 
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